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Abstract. Many data in the high energy physics are, in fact, sample means. It is 
shown that when this exact meaning of the data is taken into account and the most 
weakly bound states are removed from the hadron resonance gas, the acceptable fit 
to the whole spectra of pions, kaons and protons measured at midrapidity in central 
Pb-Pb collisions at = 2.76 TeV [Phys.Rev.Lett. 109 , 252301(2012)] is obtained. 

The invariant distributions are predicted with the help of the single-freeze-out model 
in the chemical equilibrium framework. Low pT pions and protons are reproduced 
simultaneously as well as p/n ratio. Additionally, correct predictions extend over lower 
parts of large px data. Some more general, possible implications of this approach are 
pointed out. 
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1. Introduction 

High-energy heavy-ion collisions are the tools for the creation of the deconfined phase 
(the partonic system) of the Quantum Chromodynamics (QCD) (for a wide review of the 
subject, from the theory to the experiment, see Ref. m)- The matter originated during 
such a collision, extremely dense and hot, is compressed more or less in the volume of 
the narrow disc of the ion radius at the initial moment. After then the matter rapidly 
expands due to the tremendous pressure and cools simultaneously. The evolution of the 
matter can be described in the framework of the relativistic hydrodynamics [2]. During 
expansion the matter undergoes a transition to a hadron gas phase. The hadron gas 
continues the hydrodynamical evolution, assuming that the collective behavior does not 
cease at the transition. The expansion makes the gas more and more diluted, so when 
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mean-free paths of its constituents become comparable to the size of the system one can 
not treat the gas as a collective system. This moment is called freeze-out. After then 
the gas disintegrates into freely streaming particles which can be detected. In principle, 
one can distinguish two kinds of freeze-out: a chemical freeze-out, when all inelastic 
interactions disappear and a kinetic freeze-out (at lower temperature), when also elastic 
interactions disappear. The measured hadron yields are hngerprints of corresponding 
hadron abundances present at the chemical freeze-out. The yields can be consistently 
described within the grand canonical ensemble with only two independent parameters, 
the chemical freeze-out temperature, Tch and the baryochemical potential [3]. This 
idea is the fundament of the Statistical Model (SM) of particle production in heavy-ion 
collisions. The measured pt spectra include information about the transverse expansion 
(radial flow) of the hadron gas and the temperature Tkin at the kinetic freeze-out 
|1]. However, the alternative approach to freeze-out was founded in [5l E] where the 
single freeze-out was postulated, i.e. the kinetic freeze-out coincided with the chemical 
freeze-out. This is the Single-Freeze-Out Model (SFOM). The suitably chosen freeze-out 
hypersurface and the complete inclusion of contributions from resonance decays enabled 
to correctly describe the Relativistic Heavy Ion Collider (RHIC) pt spectra. 

With the hrst data on Pb-Pb collisions at ^/snn = 2.76 TeV from CERN Large 
Hadron Collider (LHC) [TJ [8] two new problems have appeared when the SM and 
hydrodynamics were applied for the description of particle production. The predicted 
proton and antiproton abundances were much larger then measured ones |9] and low pt 
pions were underestimated laiinin]. This caused that the ratio p/tt = (p-|-p)/(7r+-|-7r ) 
was overestimated in the SM by a factor ~ 1.5 [12]. Various explanation of this ’’puzzle” 
have been invented, but all fall outside the SM. These are: (i) the incomplete list of 
resonances, there could still be undiscovered (high mass) resonances which after decays 
would increase more pion yields than proton ones, (ii) the non-equilibrium thermal 
model, with two additional parameters describing the degree of deviation from the 
equilibrium, (iii) hadronic inelastic interaction after hadronization and before chemical 
freeze-out, especially baryon annihilation, and (iv) flavor hierarchy at freeze-out, which 
could result in two different freeze-out temperatures, one for non-strange hadrons, 
another for strange hadrons (for more details and references see [T2]). 

In this work the simple generalization of the SFOM in the chemical equilibrium 
framework is postulated, where the above problems disappear naturally and whole 
results (spectra and yields) of [7] are reproduced. However, in opposite to the 
original version, all parameters of the model (thermal and geometric) are estimated 
simultaneously from the spectra. This version was successfully applied to the description 
of the hnal spectra measured at RHIC for all centrality classes in the broad range of 
collision energy [13]. The new idea introduced into the SFOM in the present work 
is to randomize one of the parameters of the model. It has turned out that the 
successful improvement is achieved only when the freeze-out temperature becomes a 
random variable and nothing is gained with the randomization of geometric parameters 
of the model. This result suggests that the temperature of the thermal system at 
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the freeze-out fluctuates significantly in the most central bin of Pb-Pb collisions at 
^/SNN = 2.76 TeV but its size remains the same. 


2. The Model 


In the SFOM the invariant distribution of the measured particles of species i has the 
form 

= j MP • “) ’ ( 1 ) 

where da^ is the normal vector on a freeze-out hypersurface, = x^/tj is the four- 
velocity of a fluid element at the freeze-out and fi is the hnal momentum distribution 
of the particle in question. The hnal distribution means that /* is the sum of primordial 
and decay contributions to the distribution. The freeze-out hypersurface is dehned by 
the equations 

Tf = H < pmax , ( 2 ) 


where the invariant time, Tf, and the transverse size, pmax, are two geometric parameters 
of the model. For the LHC energies all chemical potentials can be put equal to 
zero, so the freeze-out temperature, Tf, is the only thermal parameter of the model. 
The contribution from the weak decays concerns (anti-)protons mostly [H [H], hence 
secondary (anti-)protons from primordial and decay A(A)’s are subtracted. Fitting 
expression ([I]) to the all spectra reported in [7] (and within whole ranges presented there) 
resulted in /'f^dof=I-74: with p-value = 2 ■ 10“^^ il^dof = 235), which is unacceptable. 

Howeven the data on px spectra mn are not ’points’ but averages over all events 
in a sample 1||, the division by N^v, the number of events in the sample, means that. 
But the model prediction, Eq. ([T]), represents a quantity obtained in one collision (one 
event). Therefore, the theoretical prediction should be also an average. For that reason 
it is postulated that the expression given by Eq. o becomes a statistic (a function 
of a random variable, by dehnition it is also a random variable) and that one of the 
parameters of the model, 6 {6 = Tf,Tf or Pmax), is a random variable. Then the 
theoretical prediction is dehned as the appropriate average: 


/ dNj 

\ d^pT dy I g 


dN, 


-f{9)d9 , 


(3) 


/ d'^pT dy 

where f{9) is the probability density function (p.d.f.) of 9. This approach is more 
general but includes the standard one, if huctuations of 9 are negligible, then its p.d.f. 
is Dirac-delta like, f{9) ~ d{9 — 9o) and the average becomes the value at the optimal 
point 9o. It has turned out that only randomization of Ty improves the quality of the 
ht, randomization of Pmax or Xf does not change anything. In fact, for the technical 
reasons, not Tf is randomized but fif = I/Tf. From the statistical point of view these 


f the sample is the 0-5% centrality class here. 
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two possibilities are equivalent, because fifiTj) has a unique inverse and vice verse [T7] . 
Two p.d.f.’s are considered: log-normal 




and triangular 


1 


fWrJf,n 


T-\I3,- ft, 

r2 

0 


for \ fif — jdf |< r 
for I — /3/ |> r . 


( 4 ) 

( 5 ) 


where and a are parameters of the log-normal p.d.f. whereas fij and T are parameters 
of the triangular p.d.f., fif is the average of fif. The hrst is differentiable but has an 
inhnite tail, the second is not differentiable but has a finite range. The choice is arbitrary, 
but two general conditions should be fulhlled, a p.d.f. is dehned for a positive real 
variable and has two parameters so as the average and the variance can be determined 
independently. 

However, in both cases of p.d.f.’s, Eqs. (jl]) and (jS]), hts of expression (|3]) to the whole 
data on px spectra [7] resulted in = 1-49 with p-value = 2 • 10“® {ndof = 234), 

which is still unacceptable. 

The second assumption of the model is purely heuristic - it states that the most 
weakly bound resonances should be removed from the hadron gas. To be more precise, 
all resonances with the full width T > 250 MeV (and masses below 1600 MeV) are 
removed [U]. These are: /o(500), hi(1170), ai(1260), 7r(1300), /o(1370), 7ri(1400), 
ao(1450), p(1450), 77q(1430) and A^(1440)[^. It should be noticed that the note attached 
to /o(500) says: ’’The interpretation of this entry as a particle is controversial” [15] and 
the removal of this resonance has found the theoretical justihcation recently [T6|. The 
exclusion of only /o(500) moves fits to the boundary of the acceptance, yfi/udof ~ 1.3 ( 
p-value ~ 0.001), nevertheless according to the rigorous rules of the statistical inference 
it is still not a ’’good” £t [ 17 ]. The removed resonances are weakly bound already in 
the vacuum, with the average lifetime r < 1 fm, so they might not form in the hot 
and dense medium at all, at least in the case of central Pb-Pb collisions at extreme 
energy ^Jsnn = 2.76 TeV. Anyway, this is a heuristic hypothesis, but it works very 
well. It should be stressed at this point that both assumptions are necessary, if only the 
removal of weakly bound resonances is applied (no randomization of any parameter), the 
fit is still unacceptable, l^dof = 1-5 {p-value = 10“®). It looks like both assumptions 
(phenomena) strengthen each other. 


3. Results 

The results of hts are presented in Table [T] and depicted in Figs.[T] and [2l Predicted 
spectra are the same for positive charge particles and corresponding negative charge 

§ In fact, the hint for this assumption was the accidental observation that after up-to-date of the 
/o(500) mass to the lower one [15], the quality of fit became worse. 
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particles. This is because, when all chemical potentials equal zero, particles are 
distinguished from each other only by a mass. But particles and their antiparticles 
have the same mass and additionally patterns of decays of antiparticles are mirrors of 
patterns of decays of corresponding particles - the hadron resonance table consists of 
particles together with their antiparticles. 


Table 1. Fit results for 0-5% central Pb-Pb collisions at ^/snn = 2.76 TeV and the 
measurement at central rapidity, | j/ |< 0.5, Udof = 234. Parameters of the log-normal 
p.d.f. have no units, so their values correspond to /3f counted in MeV“^ implicitly. 


log-normal p.d.f. 

V 

Pmax 

g 

a 

E[Tf\ 

^/V[Tf] 

X^ludof 

p-value 

(fm) 

(fm) 



(MeV) 

(MeV) 


{%) 

13.80 ± 0.40 

20.48 ± 0.60 

-4.7439 ± 0.0235 

0.1764 ± 0.0090 

116.7 ± 3.0 

20.7 ± 1.6 

1.048 

29 

triangular p.d.f. 

V 

Pmax 

/?/ 

r 

E[Tf] 

^V[Tf] 

X^ludof 

p-value 

(fm) 

(fm) 

(MeV-^) 

(MeV-i) 

(MeV) 

(MeV) 


(%) 

14.42 

21.45 

0.0092482 

0.0040906 

111.6 

22.6 

1.026 

38 


Note that errors on estimates of parameters are given only for the log-normal 
p.d.f.. This is because the standard way of expressing these errors is via the Hessian 
of the fGst statistic evaluated at the estimates mi. Precisely, the inverse covariance 
matrix is given by one half of the Hessian. This requires the statistic to be at 
least twice differentiable with respect to parameters, which is not fulhlled in the case 
of the triangular p.d.f.. Both hts are acceptable and with the same quality practically. 
Predictions of both cover each other in the whole htted range and even further, they 
start to deviate from each other at pt ~ 4.5 GeV/c, see Figs. [T] and [2l As it can be 
seen from the bottom plots, predicted values stay within the two error band in general, 
only a few protons in the range 1.5-2.5 GeV/c exceed the band slightly. Additionally, 
the predictions extend over lower part of the large pt measurements [H], pions are 
reproduced up to pt ~ 4 GeV/c, kaons up to pr ~ 4.5 GeV/c and protons and 
antiprotons up to ~ 6 GeV/c. In Figs.[T] and |2] data from the large pt measurements 
are presented as one half of corresponding sums of negative and positive hadrons 
reported in Ref. [T8] . 

In Table [1] the expectation value, E[Tf], and the square root of the variance, 
\/V[Tf], of the freeze-out temperature are given. The expectation value is of the 
order of 115 MeV, whereas the square root of the variance is ~ 20 MeV. The latter 
is the measure of the fluctuations of the freeze-out temperature in the sample and is 
of the order of 20%. Since the distribution of the freeze-out temperature means the 
distribution within the sample in this work, these fluctuations reflect the changes of the 
freeze-out temperature from one event to another. Generally, these changes could be of 
the thermal and non-thermal origin. The thermal component represents fluctuations in 
an ensemble as derived e.g. in [19], i.e. it describes changes of the temperature from one 
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Figure 1. The upper panel presents spectra of positive pions, kaons and protons 
measured in central Pb-Pb collisions at -y/SNN = 2.76 TeV, data used in the fit 
are presented as error bars only [7], errors are sums of statistical and systematic 
components added in quadrature. Symbols denote one half of corresponding sums of 
negative and positive hadrons from large pt measurements m- Lines are predictions 
for the log-normal p.d.f. of /?/, dashed lines for the triangular p.d.f., in the fitted 
region they cover each other. The lower panel shows a deviation of data to the model 
in the fitted range measured in error units, (fexp — fmod)/<Jexp, where fexp(mod) is the 
experimental (model) value of the invariant yield at given px and a^xp is the error on 
fexp- In both plots circles (red online) denote pions, triangles (blue online) kaons and 
squares (green online) protons. 


system to another but prepared exactly in the same way. The non-thermal component 
expresses the variation of the freeze-out conditions event-by-event. The occurrence of 
this variation would mean that systems (events) were prepared differently in spite of the 
fact that they belong to the same centrality class. This would indicate that the 0-5% 
centrality class is inhomogeneous signihcantly. With the data available at present, it is 
impossible to distinguish between these two components. The inhomogeneity might be 
expected from the analysis of the centrality determination 121 . where the estimates of 
ranges of the impact parameter for centrality classes are given. In the 0-5% centrality 
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Figure 2. The same as Fig. [T] but for negative pions, kaons and antiprotons. 


class this parameter varies from 0 to 3.5 fm and this is the biggest spread for all classes 
considered in [5], for the rest the spread is 2 fm at most [20]. However, becanse the 
impact parameter determines the initial geometry of the collision, one conld expect 
that its variation wonld inflnence the hnal geometry of the collision, i.e. Pmax shonld 
vary from one event to another, rather. Resnlts of this work show that this is not the 
case, randomization of pmax (or Tf) has not improved the qnality of hts, they remained 
nnacceptable. Another possibility is that the variations of the hnal size of the system 
within the sample have negligible impact on the spectra. Resnlts of this model snggests 
that dnring each event a thermal system is created indeed and with approximately the 
same size at its end, however with different temperatnre. And the hnal shape of the 
spectra is the conseqnence of snmming emissions from many diherent sonrces. 

In the SFOM particle yields per nnit rapidity are given by 121 ]: 

dNi 2 

= T^Pmax'^f^i , ( 6 ) 

where n, is the thermal density of particle species i. The above expression can be 
obtained by the integration of the distribntion ([T]) over transverse momentnm. Bnt 
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here predictions are appropriate averages, so the average particle yield per unit rapidity 
reads: 

because integration over can be exchanged with the integration over (3f and p^ax and 
Tf are independent variables. Results for the yields calculated with the use of Eq. ([7]), at 
the values of parameters gathered in Table [H are given in Table |2j Excellent agreement 
with the data [7] has been achieved. Since yields are correctly reproduced, their ratios 
are correct as well. The conclusion of Ref. [16] about the enhancement of p/n ratio 
after removal of /o(500) state should be commented at this point. This is true, but 
at the same values of model parameters as before the removal, i.e. at the constant 
temperature or at the constant parameters of the temperature distribution. However 
after the removal, the expression for the model prediction, Eqs. ([T|) and ([3|), changes 
to the new one, so the f^st statistic changes as well and the new minimization is 
necessary. This new minimization estimates new values of parameters and this analysis 
proves that the new p/n ratio is smaller. 


Table 2. Midrapidity particle yields ^^||y|<o .5 ^-nd their ratios. Data are from [Tj. 


Species 

Data 

Model: 

/dm\ 

\ / Sf 

triangular p.d.f. 

log-normal p.d.f. 

7r~^ 

733 ± 54.0 

745.3 

739.2 

7r“ 

732 ± 52.0 

745.3 

739.2 

K+ 

109 ± 9.0 

106.9 

107.5 

K- 

109 ± 9.0 

106.9 

107.5 

P 

34 ± 3.0 

33.0 

32.9 

P 

33 ± 3.0 

33.0 

32.9 

Ratios 


Data 

triangular p.d.f. 

log-normal p.d.f. 

p/ir 

0.046 ± 0.003 

0.044 

0.045 

Kjr 

0.149 ± 0.010 

0.143 

0.145 


4. Conclusions 

In summary, the chemical equilibrium Single-Freeze-Out Model has been applied 
successfully to the description of the production of identihed hadrons measured at 
midrapidity in central Pb-Pb collisions at ^/snn = 2.76 TeV [7|. This has been achieved 
with the help of the more general, direct interpretation of the data and the removal of 
the most weakly bound resonances from the hadron gas. Since the chemical equilibrium 
SFOM without the above two new assumptions succeeded in the correct description of 
spectra measured at the RHIC in the broad range of collision energy 0 El DSl E3, it 
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might suggest new phenomena occurring in the most central class of Pb-Pb collisions 
at ^/SNN = 2.76 TeV. These phenomena seem to appear at the two levels: in individual 
events, where the production of identihed hadrons in each collision can be describe 
within the chemical equilibrium SFOM but with the reduced content of the hadron gas, 
and in the whole sample, causing substantial differences among collisions belonging to 
the same most central class. 

As last remarks, two more general, possible implications of the presented approach 
are pointed out. The rigorous treatment of the data ’’points” might shed some light 
on the excellent performance of the Tsallis distribution in the htting of transverse 
momentum spectra observed at the RHIC and the LHC [23]. According to Ref. [23] the 
Tsallis distribution is the mean of the exponential (Boltzmann) distribution exp {—fiE) 
weighted with a gamma distribution of fi. In the view of this work, the Tsallis 
distribution (as a mean) is the more correct equivalent of the data than the Boltzmann 
distribution. Last, but not least, a great deal of data in high energy physics are 
averages, so in any theoretical modeling (of these data) one should be aware of possible 
misinterpretations when an average is compared with a prediction for a single event. 
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